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Abstract 

A detailed discussion of the Krein's results ( applicable for solving the inverse 
scattering problem) is given with complete proofs. 

It is shown that the .S-function S(k) used in Krein's work is the S-matrix used 
in physics. 

The basic new results of the paper include the detailed description and analysis 
of an inversion algorithm based on Krein's results and a proof of its consistency, 
that is the proof that the reconstructed potential generates the same scattering 
data from which it was reconstructed. 

Numerical advantages of using Krein's method are discussed. 



1 Introduction 

The inverse scattering problem on half-axis consists of finding 

q(x) E Lii — \^q : q — q, J x\q(x)\dx < oo| 

from the knowledge of the scattering data 

S := {S(k), kj, Sj, 1 < j < J}. (1.1) 

Here 

^ ~ W 
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is the 5-matrix, f(k) is the Jost function 



f(k):=f(0,k), (1.2) 

f(x, k) is the solution to the equation 

if:=f" + k 2 f-q(x)f = (1.3) 

which is uniquely defined by the condition 

f(x,k) =e ikx + o(l), x->+oo, (1.4) 

kj > are the (only) zeros of f(x) in the region C + := {k : Imk > 0}, — k 2 are the 

■■£s + q(x) on (0, oo), Sj 



negative eigenvalues of the Dirichlet operator — 4-2 + q{x) on (0, oo), Sj > are the 



normmg constants: 

Sj = —■ , 1.5 

/('/.• ; !/'(().//.•) 

and J > is the number of the negative eigenvalues. 

For simplicity we assume that there are no bound states. This assumption is removed 
in section 4. 

This paper is a commentary to Krein's paper jl]. It contains not only a detailed 
proof of the results announced in but also a proof of the new results not mentioned in 
Q. In particular, it contains an analysis of the invertibility of the steps in the inversion 
procedure based on Krein's results, and a proof of the consistency of this procedure, that 
is, a proof of the fact that the reconstructed potential generates the scattering data from 
which it was reconstructed. Basic results are stated in Theorems 1.1 - 1.4 below. 

Consider the equation 



(I + H X )T X :=T x (t,s) + H(t-u)T x (u,s)du = H(t - s), 0<t,s<x. (1.6) 

Jo 

Equation (1.6) shows that 1^ = (I + H^H = I - (I + H x )~ l , so 

(I + H^-^I-T, (1.6') 

in operator form, and 

H={I-T x )- l -I. (1.6") 
Let us assume that H(t) is a real-valued even function 

H(-t)=H(t), H(t) G L 1 (R) fl L 2 (IR), 

/oo poo 
H(t)e ikt dt = 2 / cos(kt)H(t)dt. (1.7) 
■oo JO 
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Then (1.6) is uniquely solvable for any x > 0, and there exists a limit 

T(t,s)= \imT x (t,s):=T 00 (t,s), t,s>0, (1.8) 

x— >oo 

where T(t, s) solves the equation 

POD 

F(t,s)+ H(t-u)T(u,s)du = H(t- s), < t, s < oo. (1.9) 
io 

Given H(t), one solves (1.6), finds T 2x (s, 0), then defines 

E(x,k) - E(x,-k) 
ip{x,k):= , (1.10) 



where 

■>2x 



E(x, k) := e 



tfcx- 



1- / T 2x (s,0)e~ tks ds 
Jo 



(1.11) 



Formula (1.11) gives a one-to-one correspondence between E(x,k) and r 2a ,(s,0). 



Remark 1.1. In /i/ r 2:E (0, s) is used in place ofT 2x (s, 0) in the definition of E(x, k). By 
formula (2.21) (see section 2 below) one has r x (0, x) = T x (x,0), but F x (0,s) ^ r x (s, 0) 
in general. 

Note that 

^±ikx 



where 



and 

Note that 
where 



E(x,±k) = e ±lkx f(±k) + o(l), x^+oo, (1.12) 



roc 

f(k):=l- / T( S y ks ds, (1.13) 
./o 



T(s) := lim T x (s,0) 1=^,0), (1.14) 



#E,fc) = - Jl J . — + o(l), x-^+oo. (1.15) 



^(x, fc) = |/(fc)| sin(fcr + 5(fc)) + o(l), a: -> +00, 



The function S(k) is called the phase shift. One has S(k) = e 2i5l - k \ 
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We have changed the notations from QTJ| in order to show the physical meaning of 
the quantity (1.12): f(k) is the Jost function of the scattering theory. The function 
ip(x, k)/f(k) is the solution to the scattering problem (see equation (2.39) below). Krein 
calls 

m - f j§ die) 

the S'-function, and we will show that (1.16) is the 5*-matrix used in physics. 

Assuming no bound states, one can solve the inverse scattering problem (ISP): 

given S(k) Wk > 0, find q(x). 

A solution of the ISP based on the results of [Q consists of four steps: 

1) Given S(k), find f(k) by solving the Riemann problem (2.37). 

2) Given f(k), calculate H(t) using the formula 

Jkt •- ' 



1+H=l+ L H(t)e dt= ww (L17) 

3) Given H(t), solve (1.6) and find T x (t, s) and then r 2a; (2a;, 0), < x < oo. 

4) Define 

A(x) = 2T 2x (2x,0), (1.18) 

where 

A(0) = 2H(0), (1.18') 

and calculate the potential 

q(x) = A 2 (x) + A'(x), A(0) = 2H(0). (1.19) 

One can also calculate q(x) by the formula: 

q(x) = 2^-[T 2x (2x,0)-T 2x (0,0)]. 

Indeed, 2T 2x {2x,0) = A{x), see (1.18), £T 2x {0,0) = -T 2x {2x,0)T 2x {0,2x), see (2.22), 
and T 2x (2x,0) = T 2x (0,2x), see (2.21). 

There is an alternative (known) way, based on the Wiener-Levy theorem, to do step 

1): 
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Given S(k), find 6(h), the phase shift, then calculate the function 

2 f°° 

g(t) : = / S(k)sm(kt)dk, 

n Jo 

and finally calculate 



oo 

ikt , 



f(k) = exp / g(t)e M dk . 



o 



The potential q G L^i generates the S-matrix S(k) with which we started provided 
that the following conditions (1.20) -(1.22) hold: 



S(k) — S(-k) — S~ 1 (k), keR, (1.20) 



where 



the overbar stands for complex conjugation, 
and 



ind R S(k) = 0, (1.21) 
\\F(x)\\ Lo » (R+) + ||F(x)|| L i (R+) + \\xF'(x)\\ L i (R+) < oo, (1.22) 

1 f°° 

F(x) :=— [1 - S(k)]e ikx dk. (1.23) 

By the index (1.21) one means the increment of the argument of S(k) ( when k runs from 
— oo to +oo along the real axis) divided by 2n. The function (1.10) satisfies the equation 

ip" + k 2 ip - q{x)ip = 0, x G R+. (1.24) 

Recall that we have assumed that there are no bound states. 

In section 2 the above method is justified and the following theorems are proved: 

Theorem 1.1. Ij (1.20) -(1.22) hold, then q(x) defined by (1.19) is the unique solution 
to ISP and this q{x) has S{k) as the scattering matrix. 

Theorem 1.2. The function f{k), defined by (1.13), is the Jost function corresponding 
to potential (1.19). 

Theorem 1.3. Condition (1.7) implies that equation (1.6) is solvable for all x > and 
its solution is unique. 
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Theorem 1.4. If condition (1.7) holds, then relation (1.14) holds and T(s) := r oo (s,0) 
is the unique solution to the equation 

POO 

T(s)+ H(s-u)T(u)du = H(s), s>0. (1.25) 
Jo 

The diagram explaining the inversion method for solving ISP, based on Krein's results 
from [T], can be shown now: 

s(k) ,2 i 8 > /(*) ( W> H(t) <g r.(«, .) (tr S. ,al) r fc (!to, o) (1 4 8) a( x) (1 4?> „(*). (1.26) 

S\ S2 S3 S4 55 S6 

In this diagram s m denotes step number m. Steps S2, S4, S5 and S6 are trivial. Step si is 
almost trivial: it requires solving a Riemann problem with index zero and can be done 
analytically, in closed form. Step S3 is the basic (non-trivial) step which requires solving 
a family of Fredholm-type linear integral equations (1.6). These equations are uniquely 
solvable if assumption (1.7) holds, or if assumptions (1.20)-(1.22) hold. 

We analyze in section 2 the invertibility of the steps in diagram (1.26). Note also that, 
if one assumes (1.20)-(1.22), diagram (1.26) can be used for solving the inverse problems 
of finding q(x) from the following data: 

a) from f(k), VA; > 0, 

b) from \f(k)\ 2 , Vk > 0, or 

c) from the spectral function dp(X). 

Indeed, if (1.20)-(1.22) hold, then a) and b) are contained in diagram (1.26), and c) 
follows from the known formula (e.g., [[|, p. 256) 

{ v^ ds A > n 

* W^P' ? " (1-27) 

0, A < 0. 

Let A = k 2 . Then (still assuming (1.21)) one has: 

2k 2 1 

dp = ———dk, k>0. (1.28) 

7T \f{k)\ 2 

Note that the general case of the inverse scattering problem on the half-axis, when 

md R S(k) := v ± 0, 

can be reduced to the case v = by the procedure described in section 4 provided that 
S(k) is the S— matrix corresponding to a potential q G Li ) i(R + ). Necessary and sufficient 
conditions for such an S(k) are conditions (1.20)-(1.22) (see [4]). 
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Section 3 contains a discussion of the numerical aspects of the inversion procedure 
based on Krein's method. There are advantages in using this procedure (as compared 
with the Gelfand-Levitan procedure): integral equation (1.6), solving of which constitutes 
the basic step in the Krein inversion method, is a Fredholm convolution-type equation. 
Solving such an equation numerically leads to inversion of Toeplitz matrices, which can 
be done efficiently and with much less computer time than solving the Gelfand-Levitan 
equation (5.3). Combining Krein's and Marchenko's inversion methods yields the most 
efficient way to solve inverse scattering problems. 

Indeed, for small x equation (1.6) can be solved by iterations since the norm of the 
integral operator in (1.6) is less than 1 for sufficiently small x, say < x < xq. Thus 
q(x) can be calculated for < x < ^ by diagram (1.26). 

For x > one can solve by iterations Marchenko's equation for the kernel A(x, y): 

poo 

A(x,y) + F(x + y) +/ A(x, s)F(s + y)ds = 0, < x < y < oo, (1.29) 

J X 

where, if (1.21) holds, the known function F(x) is defined by the formula: 

1 f°° 

F(x) :=— [1 - S{k)]e ikx dk. (1.30) 

Indeed, for x > the norm of the operator in (1.29) is less than 1 (||) and it tends to 
as x — ► +oo. 

Finally let us discuss the following question: in the justification of both the Gelfand- 
Levitan and Marchenko methods, the eigenfunction expansion theorem and the Parseval 
relation play the fundamental role. In contrast, the Krein method apparently does not 
use the eigenfunction expansion theorem and the Parseval relation. However, implicitly, 
this method is also based on such relations. Namely, assumption (1.7) implies that the 
function S(k), that is, the ^-matrix corresponding to the potential (1.19), has index 0. If, 
in addition, this potential is in Li ) i(M + ), then conditions (1.20) and (1.22) are satisfied as 
well, and the eigenfunction expansion theorem and Parseval's equality hold. Necessary 
and sufficient conditions, imposed directly on the function H(t), which guarantee that 
conditions (1.20)-(1.22) hold, are not known. However, from the results of section 2 it 
follows that conditions (1.20)-(1.22) hold if and only if H(t) is such that the diagram 
(1.26) leads to a q(x) E L 1A (R + ). Alternatively, conditions (1.20)-(1.22) hold (and 
consequently, q(x) E Lx j i(M + )) if and only if condition (1.7) holds and the function f(k), 
which is uniquely defined as the solution to the Riemann problem 

$ + (&) = [1 + ^(A;)]~ 1 $_(A;), fcel, (1.31) 

by the formula 

/(*)=*+(*), (1-32) 
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generates the S'-matrix S(k) (by formula (1.16)), and this S(k) satisfies conditions (1.20)- 
(1.22). Although the above conditions are verifiable in principle, they are not satisfactory 
because they are implicit, they are not formulated in terms of structural properties of 
the function H(t) (such as smoothness, rate of decay, etc.). 

In section 2 Theorems 1.1 - 1.4 are proved. In section 3 numerical aspects of the 
inversion method based on Krein's results are discussed. In section 4 the ISP with bound 
states is discussed. In section 5 a relation between Krein's and Gelfand-Levitan's methods 
is explained. 



2 Proofs 

Proof of Theorem 1.3. If v G L 2 (0, x), then 



(v + H x v,v) = —[(v,v) L 2 {K) + (Hv,v) L 2 (R) ] (2.1) 



where the Parseval equality was used, 

PX 

v := / v (s)e iks ds, 



o 



x 

2 



v,v)= \v\ ds = (v, v) L 2 (R) . (2.2) 







Thus I + H X is a positive definite selfadjoint operator in the Hilbert space L 2 (0, x) if (1.7) 
holds. Note that, since H{t) G -^ 1 (IR), one has H(k) — > as \k\ — > oo, so (1.7) implies 

1 + H(k) > c> 0. (2.3) 

A positive definite selfadjoint operator in a Hilbert space is boundedly invertible. Theo- 
rem 1.3 is proved. □ 

Note that our argument shows that 

||(/ + ^)- 1 || L2(M) <c~\ (2.4) 

Before we prove Theorem 1.4, let us prove a simple lemma. For results of this type, see 
0- 

Lemma 2.1. The operator 

POO 

Hip := / H(t - u)(p(u)du (2.5) 
Jo 

is a bounded operator in L P (W + ), p = 1, 2, oo. 
For T x (u, s) G L 1 (R + ) one has 

poo POO 

|| / duH(t-u)T x (u, s) Iliads,) < ci( / du\T x (u, s)\) 2 . (2.6) 

J X J X 
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Proof. Let ||<p|| p := |MIlp(r+)- One has 

poo poo poo 

\\Hip\lx < sup / dt\H(t-u)\ \<p(u)\du< / = 2\\H\\ 1 \\(p\\ 1: 

«eR+ Jo Jo J-oo 

(2.7) 

where we have used the assumption H(t) = H(—t). Similarly, 

| \H(p\ |oo < 2||//||i Halloo. (2.8) 
Finally, using Parseval's equality, one gets: 



\\H<p\\l = 2n\\Hv + \\l HR) < 2ir sup \H(k)\'\\<p\\i (2.9) 

km 



where 



, + (*):={* ( t ^ °: (2.10) 

Since \H(k)\ < 2\\H\\i one gets from (2.9) the estimate: 

\\Hy\\ 2 < 2v / 2^||//||i |M| 2 . (2.11) 

To prove (2.6), one notes that 

px poo px poo 

I dt\ duH(t-u)T x (u,s)\ 2 < sup / dt\H(t - u)H(t - v)\( / |r x (u, s)|dw) 2 

JO J x u,v>x Jo J x 

poo 

< Cl ( du\T x (u,s)\) 2 . 

J x 

Estimate (2.6) is obtained. Lemma 2.1 is proved. □ 

Proof of Theorem 1.4- Define T x (t, s) = for t or s greater than x. Let w := T x (t, s) — 
T(t,s). Then (1.6) and (1.9) imply 



poo 

(I + H x )w = H{t-u)Y{u,s)du:=h x {t,s). (2.12) 

J X 

If condition (1.7) holds, then equations (1.9) and 91.25) have solutions in L X (]R + ), and, 
since sup tGR \H(t)\ < oo, it is clear that this solution belongs to L°°(IR + ) and consequently 
to L 2 (IR + ), because 1 1 1 1 2 < I M |e»| M |i- The proof of Theorem 1.3 shows that such a 
solution is unique and does exist. From (2.4) one gets 

sup Wil + H,)- 1 \\ L 2 (0<X) < c~\ (2.13) 

x>0 
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For any fixed s > one sees that sup x>y \ \h x (t, s)\ \ — > as y — > oo, where the norm here 
stands for any of the three norms L p (0, x),p = 1, 2, oo. Therefore (2.12) and (2.11) imply 



\ w \\l*(o,x) < c 2 \\ h x\\o,: 



< c 



H(t — u)T(u, s)dy 



L!(0,a:) 



H(t — u)T(u, s)dy 



L°°(0,x) 



< const ||r(u,s)|| L i (a!)0o) -> as x -> oo, 



since r(u, s) G L 1 (M+) for any fixed s > and G L 1 (M). 
Also 



s )\\l°°{0,x) ^ 2(| | /l x | l^oo^j.) + ||if x W || x,oo( 0)X )) < 

d ||r(u,s)| ||i (X;0o) + c 2 sup I |if (t - «) I || 2(M I |w| || 2(0i:r) , 



where Cj > are some constants. Finally, by (2.6), one has; 

POO 

|| w(t, s)||£2( 0jX ) < c 3 ( / |r(u, s)\du) 2 — > as x — »■ +oo. 

From (2.15) and (2.16) relation (1.14) follows. Theorem 1.4 is proved. 

Let us now prove Theorem 1.2. We need several lemmas. 
Lemma 2.2. The function (1.11) satisfies the equations 

E' — ikE — A(x)E-, E(0,k) = l, E_:=E(x,-k), 

E'_ = -ikE_ - A(x)E, E_(0, k) = 1, 

where E' = and A(x) is defined in (1.18). 
Proof. Differentiate (1.11) and get 

E' = ikE - e lkx (^2T 2x (2x, 0)e~ ik2x + 2 jf 



2x 



d(2x) 



We will check below that 



and 



dr x (t, s ) 

dx 



= -r x (t,x)r x (x, s) 



r x (t, s) = r x (x -t,x- s). 



(2.14) 



(2.15) 



(2.16) 



□ 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 
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Thus, by (2.20), 

dT 2x (s,0) 



= -T 2x (s,2x)r 2x (2x,0). (2.22) 



d(2x) 

Therefore (2.19) can be written as 

p2x 

E' = ikE - e~ ikx A(x) + A{x)e lkx / T 2x (s, 2x)e~ lks ds. (2.23) 

Jo 

By (2.21) one gets 

T 2x {s,2x) = T 2x {2x-s,0). (2.24) 

Thus 

p2x p2x 

e ikx T 2x (s,2x)e- lks ds = T 2x (2x - s,0)e lk{x - s) ds 
Jo Jo 



-ikx 



p2x 

/ T 2x (y,0)e^dy. (2.25) 
Jo 



From (2.23) and (2.25) one gets (2.17). 

Equation (2.18) can be obtained from (2.17) by changing k to —k. Lemma 2.2 is 
proved if formulas (2.20)-(2.21) are checked. 

To check (2.21), use H(—t) = H(t) and compare the equation for Y x (x — t, x — s) := ip, 

T x (x — t, x — s) + / H(x — t — u)T x (u, x — s)du = H{x — t — x + s) = H(t — s), 
Jo 

(2.26) 

with equation (1.6). Let u = x — y. Then (2.26) can be written as 

V?+ / H(t-y)<pdy = H(t-s), (2.27) 
Jo 

which is equation (1.6) for ip. Since (1.6) has at most one solution, as we have proved 
above, formula (2.21) is proved. 

To prove (2.20), differentiate (1.6) with respect to x and get: 

PX 

T' x (t,s)+ H(t-u)T' x (u,s)du = -H(t- x)T x (x,s) (2.28) 
Jo 

Set s = x in (1.6), multiply (1.6) by —F x (x,s), compare with (2.28) and use again the 
uniqueness of the solution to (1.6). This yields (2.20). 

Lemma 2.2 is proved. □ 
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Lemma 2.3. Equation (1.24) holds. 

Proof. From (1.10) and (2.17)-(2.18) one gets 

_ E" - E"_ _ (ikE - A(x)E.)' - (-ifcjL - A(x)E)' 
V ~ 2i ~ 2i ■ {Z -^ } 

Using (2.17)-(2.18) again one gets 

ij," = -k 2 ip + q(x)i), q(x) := A 2 (x) + A'(x). (2.30) 

Lemma 2.3 is proved. □ 

Proof of Theorem 1.2. The function tp defined in (1.10) solves equation (1.24) and satis- 
fies the conditions 

V>(0,Jfe) = 0, ip'(0,k) = k. (2.31) 

The first condition is obvious (in M there is a misprint: it is written that ^(0, k) = 1), 
and the second condition follows from (1.10) and (2.15): 



,.. E'(0,k) - E'_(0,k) ikE - AE_ - (ikE_ - AE) 



2i 2i 



2ik - k 
~2A ~ 



x=0 

Let f(x, k) be the Jost solution to (1.24) which is uniquely defined by the asymptotics 

f(x,k) = e ikx + o(l), x^+oo. (2.32) 

Since f(x, k) and f(x, —k) are lineraly independent, one has 

^ = Cl f(x,k) + c 2 f(x,-k), (2.33) 

where ci, c 2 are some constants independent of x but depending on k. 
From (2.31) and (2.33) one gets 

ci = /(-*), c 2 = -f(k); f(k):=f(0,k). (2.34) 

Indeed, the choice of C\ and c 2 guarantees that the first condition (2.31) is obviously 
satisfied, while the second follows from the Wronskian formula: 

f (0, k)f(-k) - /(*)/'((), k) = 2ik. (2.35) 

From (2.32), (2.33) and (2.34) one gets: 

i/>{x, k) = e ikx f{-k) + e~ lkx f{k) + o(l), x -> +oo. (2.36) 

Comparing (2.36) with (1.15) yields the conclusion of Theorem 1.2. □ 
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Invertibility of the steps of the inversion procedure 
and proof of Theorem 1.1 

Let us start with a discussion of the inversion steps 1) - 4) described in the introduction. 

Then we discuss the uniqueness of the solution to ISP and the consistency of the 
inversion method, that is, the fact that q(x), reconstructed from S(k) by steps 1) - 4), 
generates the original S(k). 

Let us go through steps 1) - 4) of the reconstruction method and prove their invert- 
ibility. The consistency of the inversion method follows from the invertibility of the steps 
of the inversion method. 

Step 1. S=>f(k). 

Assume S(k) satisfying (1.20)-(1.22) is given. Then solve the Riemann problem 



f(k)=S(-k)f(-k), 



kern. 



(2.37) 



Since ind^S^A;) = 0, one has ind^S^— k) = 0. Therefore the problem (2.37) of finding an 
analytic function f + (k) in C + := {k : ImA; > 0}, f(k) := f+(k) in C + , and and analytic 
function f-(k) := f(—k) in C_ := {k : ImA; < 0} from equation (2.37) can be solved in 
closed form. Namely, define 



[ 2m J_ 00 y-k J 



Then f(k) solves (2.37), f + (k) = f(k), /_(Jfe) = f(-k). Indeed, 

In f+(k) - In f-(k) = In S(k), keM 



(2.38) 



(2.39) 



by the known jump formula for the Cauchy integral. Integral (2.38) converges absolutely 
at infinity and In S(—y) is differentiable with respect to y, so the Cauchy integral in 
(2.38) is well defined. 

To justify the above claims, one uses the known properties of the Jost function 



roo /*oo 

/(*) = 1 + / A(0, y)e iky dy := 1 + / A(y)e iky d% 
Jo Jo 



(2.40) 



where (see 



L4(y)|<c / \q(t)\dt 



(2.41) 



9A{y) 1 /y 
dy 4 g V2 



< c 



\q{t)\dt, 



(2.42) 
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c > stand for various constants and A(y) is a real- valued function. Thus 



s( k) M Iz^ziM 1 1 of 1 



Therefore 

\nS(-k) =0 ^ 
Also 

/•oo 

f(k) = l + i A(y)ye^dy, f 
Jo 

Estimate (2.41) implies 

POO POO 

/ y\A(y)\dy<2 t\q(t)\dt 
Jo Jo 



lnS(-fc)=0(j) as \k\ -> oo, fc G 



9/ 



< oo, 



so that /'(A;) is bounded for all fc G K, as claimed. Note that 



f(-k) = f(k), keR. 

The converse step 

/(*) =► 

is trivial: 

" /(*) ■ 

If indRS 1 = then f(k) is analytic in C+, f(k) ^ in C+, f(k) = 1 + O (|) as 
fc G C + , and (2.48) holds. 

Step 2. /(Jfe) =>- if (i). 

This step is done by formula (1.17): 

The integral in (2.49) converges in L 2 {R). Indeed, it follows from (2.43) that 

2 r° / 1 



|/(A;)| 2 - 1 = -^y A'(t) sm(kt)dt + O J , |Jfe|->oo, fcel. 
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The function 

i r°° 

w(k) := - A'(t)sm(kt)dt (2.51) 
k Jo 

is continuous because tA'(t) G L 1 (IR + ) by (2.42) and w G L 2 (R) since w = o as 

\k\ — > oo, k G R. 
The converse step 

/(A;) (2.52) 

is also done by formula (1.17): Fourier inversion gives \f(k)\ 2 = f(k)f(—k), and fac- 
torization yields the unique f(k), since f(k) does not vanish in C + and tends to 1 at 
infinity. 

Step 3. H T x (s,0) T 2x (2x,0). 

This step is done by solving equation (1.6). By Theorem 1.3 equation (1.6) is uniquely 
solvable since condition (1.7) is assumed. Formula (1.17) holds and the known properties 
of the Jost function (1.4) are used: f(k) — > 1 as k —>• ±oo, f(k) ^ for k ^ 0, k G R, 
/(0) ^ since ind R S(A;) = 0. 

The converse step T x (s, 0) =>- H(t) is done by formula (1.6"). The converse step 

r 2x (2x,0) ^T x (s,0) (2.53) 

constitutes the essence of the inversion method. 
This step is done as follows: 

T 2x (2x, 0) (1 4 8) A(x) (2A7) 4 2AS) E(x, k) T x (s, 0). (2.54) 

Given A(x), system (2.17)-(2.18) is uniquely solvable for E(x,k). 

Note that the step q(x) =>- f(k) can be done by solving the uniquely 
solvable integral equation (see |§]): 

f{x, k) = e ik * + J™ Shl[k{ l~ X)] q(y)f(y, k)dy (2.55) 

with q G L 11 (R + ), and then calculating f{k) = /(0, k). 

Step 4. A(x) := T 2x (2x, 0) => q(x). 

This step is done by formula (1.19). The converse step 

q(x) =>- A(x) 

can be done by solving the Riccati problem (1.19) for A(x) given q(x) and the initial 
condition 2/7(0). Given q(x), one can find 2H(0) as follows: one finds f(x, k) by solving 
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equation (2.55), which is uniquely solvable if q G Li t i(R+), then one gets f(k) := /(0, k), 
and then calculates 2H(0) using formula (2.49) with t = 0: 

Proof of Theorem 1.1. If (1.20)-(1.22) hold, then, as has been proved in |4| (and earlier 
in a different form in |?J), there is a unique g(x) G L 11 (IR + ) which generates the given 
S-matrix S(k). 

It is not proved in that q{x) defined in (1.19) (and obtained as a final result of 
steps 1) - 4)) generates the scattering matrix S(k) with which we started the inversion. 
Let us now prove this. We have already discussed the following diagram: 

S(k) (2 4 8) f(k) (1 <4 ?) H(t) T x (s, 0) T 2x (2x, 0) (1 <4 §) A(x) q(x). (2.56) 

(1.16) 

To close this diagram and therefore establish the basic one-to-one correspondence 

S(k) q(x) (2.57) 

one needs to prove 

T 2x (2x,0) =^T x (s, 0). 
This is done by the scheme (2.54). 

Note that the step q{x) =>- A(x) requires solving Riccati equation (1.19) with the 
boundary condition A(0) = 2H(0). Existence of the solution to this problem on all of R + 
is guaranteed by the assumptions (1.20)-(1.22). The fact that these assumptions imply 
q(x) G Li ) i(IR + ) is proved in [[J and [|]]. Theorem 1.1 is proved. □ 

Uniqueness theorems for the inverse scattering problem are not given in jl| , . They 
can be found in 0-0. 

Remark 2.1. From our analysis one gets the following result: 

Proposition 2.1. If q(x) G L lil (R + ) and has no bounds states and no resonance at 
zero, then Riccati equation (1.19) with the initial condition (1.18) has the solution A(x) 
defined for all x G M + . 

3 Numerical aspects of the Krein inversion proce- 
dure. 

The main step in this procedure from the numerical viewpoint is to solve equation (1.6) 
for all x > and all < s < x, which are the parameters in equation (1.6). 
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Since equation (1.6) is an equation with the convolution kernel, its numerical solution 
involves inversion of a Toeplitz matrix, which is a well developed area of numerical 
analysis. Moreover, such an inversion requires much less computer memory and time 
than the inversion based on the Gelfand-Levitan or Marchenko methods. This is the 
main advantage of Krein's inversion method. 

This method may become even more attractive if it is combined with the Marchenko 
method. In the Marchenko method the equation to be solved is (||, [|J): 

POO 

A(x 1 y) + F(x + y)+ j A(x, s)F(s + y)ds = 0, y > x > 0, (3.1) 

J X 

where F(x) is defined in (1.23) and is known if S(k) is known, the kernel A(x, y) is to be 
found from (3.1) and if A(x,y) is found then the potential is recovered by the formula: 

= _ 2 dA{x,x) 
ax 

Equation (3.1) can be written in operator form: 

(I + F X )A=-F. (3.3) 

The operator F x is a contraction mapping in the Banach space L 1 (x, oo) for x > 0. The 
operator H x in (1.6) is a contraction mapping in L°°(0,a:) for < x < Xq, where Xq is 
chosen to that 

rx 

\H(t-u)\du<\. (3.4) 







Therefore it seems reasonable from the numerical point of view to use the following 
approach: 

1. Given S(k), calculate f(k) and H(t) as explained in Steps 1 and 2, and also F(x) by 

formula (1.23). 

2. Solve by iterations equation (1.6) for x < x < x , where xq is chosen so that the 

iteration method for solving (1.6) converges rapidly. Then find q(x) as explained 
in Step 4. 

3. Solve equation (3.1) for x > Xq by iterations. Find q(x) for x > x by formula (3.2). 

4 Discussion of the ISP when the bound states are 
present. 

If the data are (1.1) then one defines 



W = niTT if md R S(x) = -2J (4.1) 



j= i j 
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and 



W(k) = —^—w(k) if md R S(k) = -2J - 1, (4.2) 

k + ry 



where 7 > is chosen so that 7 7^ fcj, 1 < j < J. 
Then one defines 



St(k) := S{k)w 2 {k) if ind M S=-2J (4.3) 



or 

.-2 



3i(k) := iS^W (fe) if ind M S = -2J - 1. (4.4) 

Since ind]RW 2 (fc) = 2 J and ind]gW /2 (/c) = 2 J + 1, one has 

mdnS^k) = 0. (4.5) 

The theory of section 2 applies to Si(A;) and yields gi(x). From qi(x) one gets g(x) by 
adding bound states —k 2 and norming constants using the known procedure (e.g. see 

I)- 

5 Relation between Krein's and GL's methods. 

The GL (Gelfand-Levitan) method consists of the following steps (see Q, for example): 
Step 1. Given f(k), the Jost function, find 



L(x,y) :-- 



2 r°° 2 / 1 \ sin kx sin ky 

Wo \JfW~ )~k ~ 

dk (\f(k)\- 2 - 1) (cos[A;(a; - y)] - cos[k(x + y)]) 



1 

71 

:= M{x-y)- M(x + y), (5.1) 







where 



M(x):=-I dk(\f(k)\- 2 - l)cos(fcx). (5.2) 



7T 







Step 2. Solve the integral equation for K(x,y): 

K(x,y) + L(x,y)+ [ K(x, s)L(s,y)ds = 0, < y < x. (5.3) 
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Step 3. Find 



= 2 dK( f> X \ (5.4) 



Krein's function (see (1.17)) can be written as follows: 

H[t) = 2^j (\f( k )\~ 2 - 1 ) e ~ lktdk =-J o (\f(k)\- 2 -l)cos(kt)dk. (5.5) 
Thus the relation between the two methods is given by the known formula: 

M(x) = H(x). (5.6) 

In fact, the GL method deals with the inversion of the spectral foundation dp of the 
operator — + q(x) defined in L 2 (M + ) by the Dirichlet boundary condition at x = 0. 
However, if ind^S^A;) = (no bound states and no resonance at k = 0), then (see e.g. 



2k 2 dk 



dp= t -l/WI 2 ' 

' ^ 0, A<0 

so dp in this case is uniquely defined by f(k). 



A > 0, A = k 2 
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